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Abstract 



We investigate the role of the double exchange mechanism for stability of 
the metallic C-type antiferromagnetic (C-AF) phase, which was observed ex- 
"jrt ! perimentally in hole doped Lai-xSr^VOs- The double exchange model treats 

localized electrons in xy orbitals as classical 5* = 1/2 spins, which interact 

by Hund's exchange Jh with yz/zx electrons in partly filled t2g orbitals. In- 

^-^ \ eluding strong on-site Coulomb repulsion U between t2g electrons, and using 

Pj ■ slave boson method we demonstrate that C-AF and metallic phase can be 

O ' stabilized due to the specific features of the hopping in degenerate and partly 

C^, filled t2g orbitals. 

>► ! 1 Introduction 

OO , Recently undoped vanadium oxides, such as LaVOs or YVO3, have attracted much 

Cn . attention due to their puzzUng magnetic properties |l]. Large on-site Coulomb 

^T I repulsion f7 ~ 5 eV between 3d electrons leads to suppression of charge fluctuations 

vQ ■ and to the effective superexchange (SE) interaction (x J — At'^/U in the low energy 

f^ I range, where t ~ 0.2 eV is the hopping element. The cubic vanadates belong to the 

class of Mott insulators with orbital degeneracy as two d electrons, localized on the 
C^ I V'^^ ion, occupy two out of three (almost) degenerate t2g orbitals. Thus, the SE 

model includes t2g orbital degrees of freedom .2], and by itself explains the observed 
C-type antiferromagnetic (AF) phase in LaVOs, with ferromagnetic (FM) chains 
staggered in two other directions. 
Q I Even more interesting and not yet understood are the properties of doped 

O ■ Lai_a;Sra;V03, where the following phases were experimentally observed: C-AF 

and insulating one up to a; = 0.178 hole concentration, C-AF metaUic one for 
0.178 < X < 0.26, and, flnally, paramagnetic and metallic one for 0.26 < a: < 0.327 
PI- Theoretical explanation of the existence of this variety of phases has to take 
j^ ■ into account that the motion of doped holes is strongly affected by: (i) intersite SE 

interaction ex J, (m) on-site Hund's interaction oc Jh — 0.68 eV 0. Such a theo- 
retical model was recently discussed by Ishihara 0, but there the Hund's term was 
introduced in a nonsystematic way via certain projection operators, with Jh -^ cxd 
limit being assumed implicitly, and the model was solved in a way suitable for small 
doping only. Here we start with an appropriate Kondo-lattice model, systematically 



introduce the large Jh limit in which the kinetic and the Hund's parts of the Hamil- 
tonian reduce to the double exchange (DE) model |2], and solve it for any realistic 
value of doping. Since in this paper we are mainly concerned with the explanation 
of the existence of the metallic phase and the associated magnetic order, we will 
simplify the SE part of the Hamiltonian as much as possible and investigate only 
the competition of the DE and the generic SE interaction due to the occupied xy 
obitals {rixy — 1), leaving the study of the full Hamiltonian as a future problem. 
In particular we want to answer three questions: (z) can C-AF and metallic phase 
be predicted by DE, (m) if yes, then why such a phase contradicts the "intuitive" 
picture of the DE mechanism which supports homogenous FM phase, and {Hi) in 
what way the generic features of degenerate t2g orbitals help to understand the 
answers obtained in points (i) and (m), especially the differences to the case of the 
Cg electrons, playing a role in the manganites [2|. Thereby, we shall use a similar 
approximation to that employed in manganites [H] , assuming locaHzed xy electrons 
and treating them as "core spins". 

2 Double exchange Hamiltonian and the stability 
of the C-AF order 

For the sake of clarity we start with a FM Kondo-lattice Hamiltonian but we already 
include important simplifications in the SE part: 
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where: S; are "core spin" S = 1/2 operators of t2g electrons in occupied xy orbitals, 
and CTo-ct' is a vector of Pauli matrices. The fermion operators {aj^,6]^} create an 
electron with spin a in yz and zx orbitals, labeled here as a and b as these orbitals 
have no amplitude along a and b axis j2], respectively. In the relevant regime of 
large Coulomb interaction U, the fermion operators act in the restricted Hilbert 
space without double occupancies in the {a, 6} orbitals which is implemented by 
the projection operators P. Then: (z) itinerant electrons in yz or zx degenerate 
orbitals at site j can hop in the allowed [(y, z) or {x, z)] plane to the nearest neighbor 
(nn) site i, but only if there are no other electrons at site i in these orbitals (first 
and second term), {ii) Hund's exchange oc Jh aligns spins of itinerant electrons 
with "core spins" (third terms), and {Hi) occupied xy orbitals are responsible for 
AF SE interactions between nn sites in the {x,y) plane (last term), while the SE 
interactions due to itinerant electrons are not included. 

The above simplifications in the SE part oc J are motivated by the experimental 
observation that the structural transitions (which favor the orbital order and occu- 
pied xy orbitals) occur in Lai_a;Sra;V03 at lower temperature than the magnetic 
ones in the entire doping range [H]. Note that the structural transition does not 



occur in the metallic phase and then one does not expect any orbital order, similar 
to the FM metallic manganites 0. Note that the SE interactions due to the excita- 
tions of yz and zx electrons could in principle be included, and would: (?) enhance 
the AF interactions and induce their alternating orbital (AO) order in (£, ij) planes, 
and (m) would give FM interactions along the z direction, and thus further stabilize 
the C-AF phase 0. 

As in Ref. 0, we treat "core spins" S = 1/2 classically and we make a local 
unitary transformation of the basis {a|t|0), a|jO)}® {5|J0), 6|jO)} in the Hilbert 
space. The transformation makes the quantization axes of the spin of an itinerant 
electron locally parallel to the one of the "core spin", and thus it diagonalizes Hund's 
term in the spin variables. Assuming also infinite Hund's exchange Jh — > c)o, 
we leave out the itinerant electrons with antiparallel spin and what is left is a 
renormalized (a priori nonhomogenous) hopping t -^ tuij of spinless electrons. The 
absolute value of the complex parameter Uy is \uij\ — cos(%/2), with Oij being the 
relative angle between "core spins" at site i and j (the phase of Uij can be chosen 
arbitrarily at T = 0). Furthermore, we will discuss only uniform magnetic structures 
and follow the so-called uniform hopping approach (UHA) pj. One obtains then 
the following classical (i.e. Zener-type [6j) DE Hamiltonian: 

Wde = -tuy Y^ a'^^^j - t Y^ a'^fij - tu^ ^ Vfi^ - t ^ VJ), 
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where: Ux = cos{9x/2), Uy = cos{9y/2), Ox {9y) is the relative angle between spins 
in the x (y) direction and L^ is the number of sites in the crystal. The restricted 
fermion creation operators are now a\ = a]^(l — o[hi), hi = fej(l — a\ai) — a\{h\) 
creates a spinless electron at site i in yz{zx) orbital. The FM order along the z 
direction is already assumed in Eq. JSJ as it minimahzes the total energy of the 
system by maximaHzing the possibility of hopping in this direction. Note that all 
but one (spin 1/2 being classical) of the assumptions necessary to derive Eq. ^ 
from Eq. QJ seem to be reasonable approximations for cubic vanadates and for 
investigation of the problems adressed in section Q] We leave quantum treatment 
of the DE for further studies but we beheve that the main features of the DE 
mechanism in cubic vanadates can be captured by this classical Hamiltonian. 

3 Numerical results and discussion 

The ground state of Eq. (Q was found by introducing KotHar-Ruckenstein slave 
boson representation pHl- Following Ref. .JJ, we adopted this method to the orbital 
case and made the mean-field approximation (the unrestricted slave boson method, 
identical with the Gutzwiller approximation and called SBMF). It leads to the 
effective hopping terms in Hde: 

t S|Sj -^ —t c\cj , t bfbj -^ —t djdj, (3) 
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Figure 1: Magnetic phase diagram at T = obtained by solving Eq. ^ within 
SBMF approximation. Under increasing hole doping x, the C-AF order changes 
via the canted Ca** phase, with spins canted on two bonds, to a homogenous FM 
phase. 



where: q (d]^) are pseudofermion creation operators for yz (zx) orbital, and Xa {xb) 
are hole numbers in these orbitals with the constraint Xa + x^ = 1 + a;. Now Eq. 
l(2l is exactly solvable by minimizing the total energy with respect to the Ux and Uy 
variational parameters |^, which determine the magnetic structure. The resulting 
magnetic phases are shown in Fig. Q]for different values of SE J/t and hole doping 
x. We should stress that equal number of holes in both orbitals {xa = Xb) is always 
concomitant with the obtained magnetic phases, though we did not impose such 
a condition. This result follows almost entirely from the fermion properties of the 
system, namely the condition that the Fermi energies of the systems with electrons 
filling partly yz orbitals and zx orbitals should be the same. Note also that Xa = Xb 
is a neccessary condtion for the existence of the AO order. Thus, including real 
SE interactions from section Q] (z) would not contradict our results in the orbital 
sector, (ii) would only further support the stability of the C-AF phase in the spin 
sector, which is already stable for a broad range of parameters. 

Let us further analyze these results. One finds that for the broad range of 
realistic values of parameters, the DE mechanism does not win in (£, ij) planes, 
which instead have AF order and are insulating, while the holes hop only along the 
FM z direction. In this way the C-AF order characterized by the one- dimensional 
metallic behavior is stabilized. Remarkably, such a phase is stable for 0.178 < a; < 
0.26 and J/t ~ 0.178, in agreement with the experimental observations [S|. These 
results are striking since at the first sight one expects that the DE should win here 
due to the small "core spin" value S — 1/2, much smaller than in the manganites 
{S = 3/2). Instead we get results which are qualitatively similar to those obtained 
by van den Brink and Khomskii in the DE model for electron doped SrMnOa jll) . 
Though, the reason is different — this behavior here follows from very specific 
features of t2g orbitals: {i) the AF SE interactions are strictly two-dimensional 
(2D) since for the occupied xy electrons virtual hopping processes (leading to the 
SE) can occur only in the {x,y) plane — this enables metallic behavior without 
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Figure 2: Magnetic phase diagrams at T = 0, with displayed magnetic order along 
the X axis, as obtained by solving Eq. ^ and assuming: (a) no electrons in yz 
orbital, and (b) in addition no hopping along the z direction. The AF order for 
X = [marked in grey (red online)] was found by finite size scaling of the results 
obtained for large clusters. Note that the AF (FM) order is stable in the entire 
range of x along the y (z) axis [though the type of order along z direction does not 
matter in (b)]. 



any loss of the magnetic energy (hopping in the z direction does not destroy the 
magnetic order); (ii) in addition, the itinerant electrons hop between the other two 
t2g orbitals, allowing for the coexistence of the magnetic order and the metallic 
behavior, and costing no extra magnetic energy. 

To understand better this behavior we show in Fig. [21 a) stable magnetic struc- 
tures in the x direction for the systems with 2D planar (£, y) AF interactions, and 
holes being doped only into zx orbital, i.e., assuming empty yz orbitals, ^ialai = 0, 
in Eq. Q. In addition, we also studied the above case but when the hopping 
between zx orbitals occurs only along the x direction, which resembles holes doped 
into d^x^-r'^ orbitals in the manganites [stable phases are displayed in Fig. El^b)]. 
One finds that for this Cg-type orbital the canted (Ca) structure is stabilized in a 
broad range of parameters, in contrast with t2g case where AF order in the plane 
can be conserved. This is one of the reasons why the AF order along the z direction 
can be so easily destabiHzed by doping in LaMnOa, leading to the FM metallic 
phase |H]- By looking at Fig. El^a), one can also realize better the role of orbital 
degeneracy: even if we had a system with planar AF interactions and we doped it 
with holes into, e.g., zx orbitals, we would not get the AF plane stable for a realis- 
tic range of parameters! Instead, one AF bond has to be broken. This tendency is 
suppressed if we have two degenerate orbitals as in Eq. ^ . 

In summary, we have shown that the existence of the metallic C-AF phase in 
doped Lai_a;Sra;V03 can be explained by the DE model with extremely simplified 
SE interactions. The stability of this phase is due to the specific features of the 
t2g orbitals — namely due to strictly 2D and strictly fiavor conserving hopping 
between these degenerate orbitals J2l- In our view these are very prospective results 
and, in connection with those of Ref. j^], they suggest that including the real 



SE interactions could indeed explain the observed insulating phase and the metal- 
insulator transition which occurs at finite doping. 
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